The Chemical Enrichment History of the Fornax Dwarf 
Spheroidal Galaxy from the Infrared Calcium Triplet 1 
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Abstract 

Near infrared spectra were obtained for 117 red 
giants in the Fornax dwarf spheroidal galaxy with 
the FORS1 spectrograph on the VLT, in order to 
study the metallicity distribution of the stars and 
to lift the age-metallicity degeneracy of the red gi- 
ant branch (RGB) in the color-magnitude diagram 
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(CMD). Metallicities are derived from the equiv- 
alent widths of the infrared Calcium triplet lines 
at 8498, 8542, and 8662 A, calibrated with data 
from globular clusters, the open cluster M67 and 
the LMC. For a substantial portion of the sample, 
the strength of the Calcium triplet is unexpect- 
edly high, clearly indicating that the main stellar 
population of Fornax is significantly more metal- 
rich than could be inferred from the position of 
its RGB in the CMD. We show that the relative 
narrowness of the RGB in Fornax is caused by the 
superposition of stars of very different ages and 
metallicities. 

The region of parameter space occupied by 
the most metal-rich red giants in Fornax -young, 
metal-rich and luminous- is not covered by any 
of the calibrating clusters, which introduces un- 
certainty in the high end of the metallicity scale. 
Using published theoretical calculations of Call 
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triplet equivalent widths, we have investigated 
their sensitivity to luminosity, age, and metallic- 
ity. The correlation between absolute I magnitude 
and Call strength appears to be only slightly af- 
fected by age variations, and we have used it to 
estimate the metallicities of the Fornax stars. 

The metallicity distribution in Fornax is cen- 
tered at [Fe/H]=— 0.9, with a metal-poor tail ex- 
tending to [Fe/H]~ —2. While the distribution to 
higher metallicities is less well determined by our 
observations, the comparison with LMC data in- 
dicates that it extends to [Fc/H]~ —0.4. By com- 
paring the metallicities of the stars with their po- 
sitions in the CMD, we have derived the complex 
age-metallicity relation of Fornax. In the first few 
Gyr, the metal abundance rose to [Fe/H]~ —1.0 
dcx. The enrichment accelerated in the past ~ 1-4 
Gyr to reach [Fe/H]~ -0.4 dcx. More than half 
the sample is constituted of star younger than ~ 4 
Gyr, thus indicating sustained recent star forma- 
tion in Fornax. These results are briefly compared 
to the theoretical predictions on the evolution of 
dwarf galaxies. They indicate that the capacity of 
dwarf spheroidal galaxies to retain the heavy ele- 
ments that they produce is larger than expected. 

1. Introduction 

During the last decade, the deep, wide field 
color-magnitude diagram (CMD) studies and 
spectroscopy of the dwarf spheroidal (dSph) satel- 
lites of the Milky Way have beautifully confirmed 
early hints on the fact that these galaxies had 
experienced more than one episode of star for- 
mation and nucleosynthesis (Zinn 1981, and ref- 
erences therein). Indeed, we find almost every 
imaginable evolutionary history in this sample of 
galaxies, from the extreme case of Leo I, which 
has formed over 80% of its stars in the second half 
of the life of the Universe (Gallart et al. 1999), to 
intermediate cases like Carina (Smcckcr-Hane et 
al. 1996; Hurley-Keller, Mateo & Nemec 1998) 
and Fornax (Stetson, Hesser & Smecker-Hanc 
1998; Buonanno et al. 1999), with prominent 
intermediate-age populations, to predominantly 
old (~ 10 Gyr old) systems like Sculptor (Hurley- 
Keller, Mateo & Grebel 1999), Draco (Aparicio, 
Carrera & Martinez-Delgado 2001), Ursa Minor 
(Mighell & Burke 1999; Carrera et al. 2002) and 
Leo II (Mighell & Rich 1996). 



The very extended star formation history 
(SFH) in several of these galaxies seems to be 
independent of their total mass: Carina and 
Leo I are among the least massive dSph in the 
Milky Way system, with virial masses around 
2 x 1O 7 M0, while Fornax is the most massive 
one, with 7 x 1O 7 M0, see Mateo (1998). Their 
metal content, however, is directly related to their 
total luminosity, and presumably, to their total 
mass (Caldwell et al. 1998 and references therein): 
Leo I and Carina have low metallicities (Gallart 
et al. 1999; Smecker-Hane et al. 1999), while 
Fornax seems to have a relatively high metallic- 
ity and a large metallicity dispersion (Saviane et 
al. 2000; Tolstoy et al. 2001; and much more 
prominently in this paper). The metal content, 
therefore, seems to be not as much related to the 
SFH than to the ability of these systems to retain 
the produced metals, which may have to do with 
the effect of SNe on the interstellar medium and 
the depths of their potential wells (e.g., Mac-Low 
& Ferrara 1999; Ferrara & Tolstoy 2000). 

The nearest galaxies offer an excellent opportu- 
nity to test these correlations and their theoretical 
interpretation, since their SFH can be derived with 
great accuracy from deep CMDs, their masses can 
be calculated from the velocity dispersion of their 
stars, and the metallicities and metallicity distri- 
bution can be obtained from spectroscopy of their 
individual stars. It is important to emphasize the 
need to replace the common technique of infer- 
ring the mean and the dispersion in metal abun- 
dance from the mean color and color spread of the 
RGB in the CMD under the assumption that the 
stars in these galaxies are very old. Because of the 
age-metallicity degeneracy in the position on the 
RGB, this method can lead to significant errors 
if the galaxies contain, as some do, sizable popu- 
lations of intermediate-aged stars. Spectroscopic 
metallicity determinations exist for a small sam- 
ple of the nearest dSph using the Ca II triplet 
and other metallicity indicators in low dispersion 
spectra (Draco: Lchnert et al. 1992, and refer- 
ences therein; Winnick & Zinn, in prep.; Sextans: 
Da Costa et al. 1991; Suntzeff et al. 1993; Carina: 
Smecker-Hane et al. 1999; Sculptor: Tolstoy et al. 
2001; Fornax: Tolstoy et al. 2001). All of the stud- 
ies that have observed more than a few stars have 
discovered that substantial metallicity dispersions 
are present, even though the mean metallicity of 
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each of these low-mass systems is low. High dis- 
persion abundance determinations of a number of 
elements require necessarily an 8m-class telescope, 
and exist for a few stars in the nearest north- 
ern dSph, namely Draco, UMi and Sextans, for 
which Keck HIRES spectroscopy has been pos- 
sible (Shetrone, Cote & Sargent 2001; Shetrone, 
Bolte & Stetson 1998) and in the nearest southern 
dSph, Sculptor, Fornax, Carina and Leo I, using 
UVES at the VLT (Shetrone et al. 2003; Tol- 
stoy ct al. 2003). These studies have measured 
just a few stars in each galaxy, and confirm the 
low resolution spectroscopy results in the cases 
where there is overlap, in terms of mean metal- 
licity and existence of a dispersion in metallicity. 
Among the new information that they provide is 
the abundance patterns among the dSph galaxies, 
which are quite uniform, indicating similar nucle- 
osynthetic histories and presumably similar IMF. 
The [a/Fe] element ratios are however in general 
lower in the dSph than in the Galactic Halo, for 
the same range of metallicity. 

In this paper we present spectroscopy and metal 
abundances for a large number of stars in the For- 
nax dSph galaxy. Fornax is a particularly inter- 
esting galaxy in the context of the study of chemi- 
cal enrichment processes in dSph galaxies, since it 
seems to be among the few Milky Way dSph satel- 
lites (the other outstanding one being Sagittar- 
ius, see Layden & Sarajedini 2000 and references 
therein) that has been able to retain substantial 
amounts of metals during its evolution, as hinted 
by the width of its RGB (Saviane ct al. 2000) 
and confirmed and shown to be even more ex- 
treme by the spectroscopic work presented in this 
paper. Tamura, Hirashita & Takeuchi (2001) con- 
sider it as the only galaxy in our immediate neigh- 
borhood with properties comparable to the more 
massive dwarf elliptical (dE) Andromeda compan- 
ions. Fornax is indeed one of the most massive of 
the Milky Way satellite dSph galaxies, with total 
mass 6.8 x 10 7 M© (Mateo 1998), as inferred from 
the central velocity dispersion of its stars (Mateo 
ct al. 1991), and it is one of the two (again with 
Sagittarius) containing its own system of globular 
clusters. The existing deep CMDs (Stetson et al. 
1998; Buonanno et al. 1999; Saviane et al. 2000) 
show a blue horizontal-branch and a well popu- 
lated red-clump indicating a substantial popula- 
tion of old and intermediate-age stars respectively, 



and a relatively bright main sequence which must 
contain stars as young as a few hundred million 
years. Particularly puzzling is the fact that no 
HI gas is found in deep VLA observations near 
the center of the galaxy (Young 1999), where the 
youngest stars are found. The current data cannot 
definitely exclude, however, the presence of HI in 
the outer parts of Fornax, where it may reside if 
it was ejected by the last event of star formation. 

To understand in more detail the evolution 
of Fornax, one would like to measure its age- 
metallicity relation (AMR) and its SFH. This can 
be done if the age-metallicity degeneracy of the 
RGB and other features of the CMD can be bro- 
ken by observations that yield the abundances of 
the stars independently of their ages. The infrared 
Ca II triplet has been shown to be a useful metal- 
licity indicator for metal-poor stars (see e. g. Rut- 
ledge et al. 1997a, hereafter RHS, and references 
therein). While most previous investigations that 
have used the infrared Ca II triplet as a metal- 
licity indicator have worked on globular clusters 
or other very old star systems, a few have col- 
lected observations of red giants in significantly 
younger stellar populations and have shown that 
the metallicities obtained in this way are consis- 
tent with other estimates (Olszewski et al. 1991; 
Suntzeff et al. 1992; Da Costa & Hatzidimitriou 
1998). This previous work motivated us to secure 
spectra at the Ca II lines of a sample of 117 stars 
in Fornax that lie near the tip of its RGB. We sus- 
pected that this sample would contain stars cov- 
ering nearly the entire ranges in age and in metal- 
licity. To demonstrate that these observations are 
not compromised by an age-metallicity degener- 
acy, we discuss the sensitivity of Ca II strength to 
age and metallicity variations using the theoreti- 
cal calculations by Jorgensen, Carlsson & Johnson 
(1992), and we show that the observations of star 
clusters are consistent with the theoretical expec- 
tations. We then derive the metallicities of the 
Fornax stars from a new calibration of the Ca II 
triplet, estimate the AMR and discuss the SFH of 
this galaxy. 

2. Observations and Reductions 

The observations were obtained on the night 
of 1st December 1999, using the FORS1 instru- 
ment installed on the Cassegrain focus of the 
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UT1 telescope (first 8.2-m unit of the ESO Very 
Large Telescope) at Paranal, Chile. FORS1 is 
a low-dispersion spectrograph that allows the si- 
multaneous spectroscopy of up to 19 objects in 
a 6.8'x6.8' field of view through a system of ad- 
justable slitlets. Grism GRIS-600I+15 was used 
with order separation filter OG590+72, resulting 
in a central wavelength of 7940 A, a dispersion of 
1.06 A per pixel (44 A/mm), and a resolution of 
~ 1530. The spectral range covered was in general 
7000-9000 A, with variations depending on the ex- 
act position of the slitlct in the field. In position- 
ing the slitlets, we ensured that the whole spectral 
region of interest to measure the Ca II triplet fea- 
ture was always covered. Calibration exposures 
were taken for each setup in the morning follow- 
ing the observations with a Ne-Ar lamp. The ex- 
cellent seeing conditions during the observations 
-most of the time below 0.8" seeing- allowed us 
to use slitlet widths of 0.7 arcsec. 

Seven fields were observed in Fornax, with an 
average of 17 targets per field. For each field, two 
20-minute exposures were acquired on two posi- 
tions offset vertically by about 3 arcsec (~10 pix- 
els). This procedure allows the direct subtraction 
of the sky from the unextracted spectra, in spec- 
troscopic analogy to the method widely used in 
infrared photometry. The method is a simplified 
version of the "va-et- vient" procedure presented 
by Cuillandre et al. (1994). It is much more ef- 
fective at removing the night sky emission lines 
than the conventional method of interpolating be- 
tween sky spectra that are recorded on either side 
of the object spectrum. Short exposures were also 
acquired on selected red giants of three nearby 
globular clusters, NGC 2298, NGC 3201 and NGC 
7099=M 30, for the purpose of the metallicity cal- 
ibration. A single exposure with an exposure time 
of 2-5 minutes was obtained for each cluster. Be- 
cause red giants near the tip of the RGB are much 
less dense in the plane of the sky in the globular 
clusters than in Fornax, only a small fraction of 
the slitlets could be used in the cluster exposures. 

2.1. Target selection 

The program targets were selected in the up- 
per part of the RGB in Fornax within the central 
11 arcmin of the galaxy. We used the B and R 
band photometry of Stetson ct al. (1998) for this 
selection. For the purpose of the target selection 



only, the / magnitude was extrapolated from the 
B and R using a relation fitting globular cluster 
data: I ~ R - \{B - R) + 0.35. The V and I 
values used later in the article and listed in Ta- 
ble 2 are from the new photometry of Fornax by 
Gallart ct al. (in prep.). The targets were se- 
lected with linear limits on color and magnitude 
encompassing all the upper part of the Fornax 
RGB down to I = 18.35: 16.60 < I < 18.35 and 
-3.24(5-5) + 23.02 < I < -2.61(5 - R) + 23.52. 
The / magnitude was used to define the selection 
because it corresponds to the wavelengths of the 
Ca II triplet region. The color selection limits 
were designed in order to include the whole width 
of the Fornax RGB, including potential outliers 
on the blue side, while excluding regions obviously 
dominated by field stars. 

This selection produced 2534 potential targets. 
The actual targets were picked according to their 
position in the field of the instrument, in order 
to fill optimally the 19 slitlets available for multi- 
object spectroscopy on FORS1. The position of all 
observed targets is shown on Fig. 1 on the plane 
of the sky and Fig. 2 in the CMD. 

While age, metallicity or kinematic criteria 
were not used in the selection, some age and metal- 
licity selection is however implicit in the color and 
magnitude limits. The upper part of the RGB 
is not a totally representative sample of a given 
stellar population. In Section 3.2 we attempt to 
infer the metallicity distribution of the complete 
population from the distribution of our sample. 

2.2. Spectra reduction 

The spectra were reduced using the noao pack- 
ages twodspec.apextract and onedpec in IRAF 5 . 
The procedure was somewhat different for the pro- 
gram stars, for which two offset frames were ob- 
tained, and for the calibrating globular clusters, 
for which only one frame of short exposure was 
obtained. 

Fornax frames 

The frames were flat-fielded with lamp flatfields 
taken after the night with the same spectrograph 

5 IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Uni- 
versities for Research in Astronomy, Inc., under contract to 
the National Science Foundation. 
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and slitlet configuration than the program frames. 
A variable width at the side of each slitlet was re- 
moved to avoid reflections. In three cases, namely 
targets 217, 302 and 714 6 , this procedure removes 
part of the spectra. 

Cosmic rays were removed with the IRAF task 
cosmicrays, and the two offset frames were sub- 
tracted to obtain an image with a positive and 
negative spectrum, some 10 pixel apart, for each 
slitlet. The positive and negative spectra were ex- 
tracted with a width of 8 pixels and calibrated in 
wavelength using Ne-Ar lamp exposures taken at 
the end of the night and the spectral line wave- 
lengths listed in the FORS user manual (In the 
case of frame 5, the wavelength calibration expo- 
sure was corrupted and re-done several weeks after 
the observation). The residual wavelength scatter 
was smaller than 0.1 A. Attention was given to 
the stability of the calibration between 8400 and 
8700A -the region of the Ca II triplet lines- allow- 
ing for possible non-linear distortions on the edges 
of the spectral range. The two offset, wavelength 
calibrated spectra, one positive and one negative, 
were then subtracted to produce the final spec- 
trum. 

This procedure allows a very efficient removal 
of the sky background. Most of the sky cancels 
out in the direct subtraction of the two images, 
independently of the variations in pixel response. 
It is therefore unaffected by the errors in the flat- 
fielding. The only sky light left in the subtracted 
images is that due to variation of its intensity be- 
tween the two 20-minutes exposures, which is re- 
moved when the two extracted, wavelength cali- 
brated spectra are subtracted again. An example 
of the raw, sky and final object spectra is given in 
Fig. 3. 

The reduction procedure was checked by ex- 
tracting the spectra of all the objects of frame 4 by 
two authors independently, using slightly different 
parameters. The resulting spectra were essentially 
identical. Varying the extraction width between 1 
and 10 pixels was also tested on frame 4, and 8 
pixels was judged to be a reasonable value to ap- 
ply to all spectra. While marginally better results 

'We use the following nomenclature for our Fornax targets: 
the hundreds digit indicates the FORS frame, and the fol- 
lowing digits number the slitlets in the frame, from top to 
bottom on the CCD. The coordinates of the targets are 
given in Table 2. 



could have been obtained by tailoring the extrac- 
tion width to each spectrum, this would have in- 
troduced a subjective element regarding what con- 
stituted a "bad feature" . Consequently, we used 
the same extraction width for all spectra. 

Globular cluster calibration frames 

Only a single exposure was obtained of the glob- 
ular cluster red giants that serve to calibrate the 
pseudo-equivalent widths, because they are much 
brighter than the sky background. Their spectra 
were extracted after de-biasing and flat-fielding. 
The sky was removed by a linear fit to the back- 
ground perpendicular to the direction of the dis- 
persion. 

We observed several RGB stars in three globu- 
lar clusters, NGC 2298, NGC 3201 and NGC 7099, 
which were selected from the photometry of Al- 
caino & Lillicr (1986), Lee (1977) and Dickens 
(1972) respectively. We chose preferentially stars 
whose Ca II triplet equivalent widths were mea- 
sured by Rutledge et al. (1997a, hereafter RHS). 
The exposure times for the clusters were 300 sec, 
180 sec and 120 sec respectively. Table 1 lists the 
information on the measured stars. The Star ID 
and the B, V magnitudes are taken from the above 
references. The equivalent widths of the three 
Ca II lines and their sigma are listed in columns 5- 
10. The combined equivalent widths as measured 
in this work and its sigma are listed in columns 11- 
12, and the same information as measured by RHS 
is listed in columns 13-14. Stars DP91, DP23, and 
DP24 in NGC 7099 were not observed by RHS, 
and the values listed are the result of transform- 
ing the observations of Suntzeff et al. (1993) to 
the RHS scale using the equation given by RHS. 

2.3. Radial velocities 

Radial velocities were obtained by fitting each 
of the three Ca II lines with a Gaussian profile 
in IRAF. Velocities were then computed using the 
measured line centers and the known lab wave- 
lengths (taken from RHS). Whenever possible, all 
three lines were used to compute the final averaged 
velocity for an individual star. 

The low resolution of the spectra limits the 
usefulness of the radial velocity to the identifica- 
tion of clear outliers as likely non-members of For- 
nax, without allowing a more detailed kinematical 
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study. Objects 109, 214 and 418 have discordant 
velocities and are suspected non-members. The 
velocity of 216 is also discordant, but the spec- 
trum is too noisy for the indication to be reliable. 
In the following analysis, the first three objects 
are excluded, while the fourth falls much below 
our spectral quality criteria and is excluded any- 
way (Par. 2.4). 

2.4. Ca II equivalent width 

Pseudo equivalent widths were calculated from 
the spectra using passbands in and around the 
Ca II triplet lines. To place our measurements 
on the system of RHS, we followed closely their 
prescription for measuring the pseudo-equivalent 
widths. We used the same functional form, the 
Moffat function of exponent 2.5, for fitting the 
lines and the same wavelength intervals for con- 
tinuum passbands. For the measurement of the 
combined equivalent width we used their defini- 
tion of the weighted sum of the widths, £Ca, which 
they defined as £Ca= 0.5u>l + w2 + 0.6w3, where 
wl,w2 and w3 are the equivalent widths of indi- 
vidual Calcium lines (A 8498, 8542, 8662A, respec- 
tively). We observed 19 red giants in the globular 
clusters NGC 2298, NGC 3201, and NGC 7099 for 
which values of ECa were measured by RHS or 
could be transformed to their system (see Table 
1). The mean difference in the sense ours minus 
theirs between the values of SCa is 0.04, with a 
standard deviation of 0.11. Since these deviations 
are small compared to the average measuring er- 
rors (0.16 and 0.10 for our and RHS values, respec- 
tively) , no corrections appear necessary to put our 
measurements on the RHS system. 

As a consistency check, the equivalent width 
were also calculated by fitting a Gaussian profile 
on the spectra with the " splof task in IRAF, af- 
ter the spectra continuum were flattened by di- 
vision with a 6-th order polynomial. The Gaus- 
sian and Moffatian values were found to be in 
excellent agreement. A regression between the 
equivalent widths calculated with the two methods 
gives ECa Mo // = 1.005SCa Go « ss +0.014A, er=0.12 
A. The passband value was adopted except when 
there was a visible spike in the continuum near the 
Calcium lines. In that case the Gaussian fit was 
preferred. In the case of frame 5, only the Gaus- 
sian values were computed. Table 2 indicates when 
Gaussian fitting was chosen. 



The signal-to-noise ratio per pixel was evalu- 
ated by computing the dispersion of the normal- 
ized background in three different wavelength in- 
tervals around the Ca II lines : A = 8390-8489 A, 
8555 - 8650 A, 8695 - 8794 A. Spectra were in- 
cluded in the analysis when the signal-to-noise ra- 
tio in these intervals was higher than 20 (mctallic- 
ities are nevertheless given in Table 2 down to a 
signal-to- noise ratio of 14). Note that this criteria, 
which takes into account the noise on the side of 
the Ca II lines only, will not produce any bias at 
selectively eliminating spectra of lower metallicity 
stars. Stars rejected by this criteria or for other 
reasons mentioned in the notes are shown as open 
dots in Figure 2 (23 out of 118 stars). In Figure 4, 
three sample spectra are shown. The final uncer- 
tainty of the combined Ca II equivalent widths is 
in the range 3-6% for our objects. 

3. Calibration and results 

3.1. Metallicity calibration of the Ca II 
triplet 

As discussed in the introduction, previous stud- 
ies have shown that the strengths of the infrared 
Ca II triplet lines can be used as a metallicity in- 
dicator for old, metal-poor stars, such as the ones 
in globular clusters. It is less clear that these lines 
provide an accurate metallicity ranking for red gi- 
ants of near solar composition and for ones span- 
ning a wide range in age. Since Fornax contains 
red giants with these properties in addition to very 
old, metal-poor ones, it is imperative that wc in- 
vestigate the behavior of the Ca II lines in some 
detail. 

In the Appendix, we use theoretical calcula- 
tions of Ca II line strength as a function of ef- 
fective temperature and log g, together with theo- 
retical isochrones, to investigate the sensitivity of 
the Call lines to age and metallicity 

Most previous investigations that used SCa as 
a metallicity indicator have plotted it against V- 
V(HB) because this quantity is independent of the 
reddening and the distance modulus of the sys- 
tem. This quantity works fine for very old stel- 
lar populations that have well-defined horizontal 
branches and small or virtually non-existent in- 
ternal metallicity spreads, such as most globular 
clusters. The metallicity calibration of the plot of 
SCa vs. V-V(HB) corrects implicitly for the sys- 
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tematic variation in HB luminosity with metallic- 
ity. It is much less clear that V-V(HB) is a useful 
parameter for a galaxy like Fornax that possesses 
wide ranges in age and metallicity. In the case 
of an intermediate-age star cluster, one can make 
the necessary correction to V-V(HB) because the 
ages of the red giants are known (e.g. Da Costa 
& Hatzidimitriou 1998). But in Fornax, because 
the age of an individual red giant cannot be deter- 
mined, one does not know if its V-V(HB) should 
be corrected for age. 

Since the distance modulus of Fornax is well 
determined and differential reddening is negligi- 
ble, it is possible to determine unambiguously 
the absolute magnitude of every star. If plots of 
SCa against absolute magnitude yield a metallicity 
ranking that is reasonably insensitive to age, then 
it docs not matter that the ages of the red giants 
span a large range and that the ages of individual 
stars cannot be determined. 

Fig. 5 shows the Fornax and cluster data in the 
projections relevant to the metallicity calibration. 
The discussion in the Appendix indicates that: 

i) When clusters of different ages are consid- 
ered, indicators primarily sensitive to gravity, like 
absolute magnitude, are expected to be much bet- 
ter for the SCa calibration than indicators sensi- 
tive primarily to temperature, for example (V-I). 
The relatively young and metal-rich open cluster 
Mil 7 provides an important check on the age sen- 
sitivity. Suntzeff et al. (1992, 1993) measured the 
infrared Ca II lines in the spectra of 21 red gi- 
ants in Mil, which we have transformed to SCa 
using the transformation equations in RHS. Mil 
can be compared to the older, similarly metal-rich 
calibrating cluster M67 8 . Fig. 5 confirms the im- 
portance of gravity effects on the Ca II triplet 
equivalent width. Mil and M67 occupy widely 
different positions in the (V-I) vs SCa plot and 
the Mi vs (V-I) plot, even though their metallic- 
ity is similar, because their age difference causes 
large gravity differences at a given temperature. 

ii) An abundance calibration of the Ca II triplet 
using the Mj magnitude is somewhat less sensitive 
to age than one using My- Our calibration of Wo 
(defined below as SCa at M/=0) with [Fe/H] on 

7 Age ~ 0.25 Gyr (Sung ct al. 1999), metallicity 

[Fc/H]=+0.10 ± 0.04 (Gonzalez & Wallerstcin 2000) 
8 Age ~ 4 Gyr, metallicity [Fe/H]~ (Richer et al. 1998). 



the Carretta & Gratton scale (see below) yields a 
mean [Fc/H] of +0.02 ± 0.03 for Mil, which sug- 
gests a remarkably low sensitivity of the SCa - Mj 
diagram to age. We have also derived a calibration 
using My instead of Mj. This calibration yields a 
mean [Fe/H] of -0.21 ± 0.03 for Mil. 

For these reasons, in the following analysis, we 
derive a calibration based on Mj because it should 
produce more accurate results for the youngest 
stars in Fornax. It is very important to empha- 
size, however, that the fairly good results obtained 
for Mil and M67 do not guarantee that the age 
effect is negligible in the case of the most metal- 
rich stars in Fornax. These stars lie in a region of 
parameter space that is not covered by the cali- 
brators (i. e. bluer than M67, brighter than Mil, 
and younger than globular clusters). In partic- 
ular, the models tend to predict a curvature of 
the isometallicity lines in the Mj vs. W plane for 
the brightest and most metal-rich red giants (see 
Fig. 16 in the Appendix). This curvature is not ob- 
served in the relatively low-luminosity M67 data. 
If present at higher luminosity, it would lead to an 
age-dependent shift of the derived metallicities for 
the metal-rich end of the Fornax data. It would 
also, however, cause a systematic trend in the re- 
covered metallicities as a function of luminosity, a 
trend that is not observed in our data (see Fig. 6). 

In the following paragraphs, we first present a 
metallicity calibration in the SCa - Mr plane based 
on the calibrating clusters, and then discuss how it 
should be modified for the youngest Fornax stars. 
Our cluster calibration of the SCa - Mr relation is 
based on the observations of Ca II line strengths in 
11 globular clusters and M67 from Olszewski et al. 
(1991), Armandroff & Da Costa (1991), Suntzeff et 
al. (1992, 1993), Da Costa & Armandroff (1995) 
and RHS. When it was necessary, the observed 
values of W(Ca II) were transformed to SCa us- 
ing the equations in RHS. The I photometry of 
the stars was taken from Da Costa & Arman- 
droff (1990), Ortolani et al. (1990, 1992), Ortolani 
(priv. comm.), Alcaino & Liller (1986), Alcaino, 
Liller & Alvarado (1989), Lloyd-Evans (1983) and 
Janes & Smith (1984). For a minority of the stars, 
I-band photometry was not available in the liter- 
ature, and it was necessary to derive approximate 
values from transformations between (B-V)o and 
(V-I)o- This was accomplished using stars that 
had been measured in both colors in the clus- 
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ters. For the very metal-poor cluster NGC7099, 
the (B-V)o vs (V-I)o transformation derived by 
Zinn & Barnes (1996) was used. The uncertain- 
ties introduced by these transformations are small 
compared to the other sources of error in the cali- 
bration. The reddenings and V magnitudes of the 
horizontal branches of the globular clusters were 
adopted from either the photometric catalogs or 
from the 1999- version of the compilation by Harris 
(1996). The distance moduli of the globular clus- 
ters were calculated from the V(HB) values using 
the My(HB) relation calculated by Demarque et 
al. (2000), which takes into account the HB mor- 
phologies of the clusters as well as their [Fe/H]. 
The apparent distance modulus and reddening of 
M67 were taken from Twarog ct al. (1997). 

While the ECa data for all 12 clusters are well 
correlated with Mj, the correlation for three clus- 
ters are clearly inferior to the other nine. For one 
of the three clusters (NGC1851), ECa was avail- 
able for only 5 stars that span a relatively small 
range of Mj. The other two clusters (NGC6528 
& NGC6553) lie in very crowded and heavily red- 
dened fields, which may have produced larger than 
average errors in line strength and Mi. The tight 
sequences produced by the other nine clusters are 
well approximated by straight lines of nearly the 
same slope, and there is no firm evidence of sys- 
tematic change in slope with [Fe/H]. The aver- 
age slope is —0.48 ± 0.02, which is also not far 
from the values obtained with data for NGC1851, 
NGC6528 and NGC6553. Following previous in- 
vestigators, wc define a reduced equivalent width 
(Wo) to remove the increase in ECa with increas- 
ing luminosity. 

W a = E Ca + 0.48 Mj 

Note that this Wo is ECa at Mj=0, and therefore 
not the same as W at V— V(HB)=0 as defined by 
other authors. 

To calibrate Wq in terms of [Fe/H], we adopted 
for eight of the globular clusters the value of [Fe/H] 
that Carretta & Gratton (1997, CG) and Carretta 
et al. (2001) measured from high dispersion spec- 
tra of red giants. Carretta et al. (2001) consider 
their measurements of the most metal-rich clusters 
NGC 6528 and 6553 to be on the same scale as the 
more metal-poor ones that were measured by CG. 
For NGC1851 and M2, which were not measured 
by these teams, we selected the values that RHS 



obtained from a calibration, on the CG scale, of 
the plot of V(HB)-V against ECa. Several stud- 
ies have shown that M67 has a solar composition 
within the errors, and we adopted [Fe/H]=0.00 ± 
0.10 from Twarog et al. (1997). Our calibration 
of Wo with [Fc/H] is shown in Figure 7. The solid 
curve is the expression: 

[Fe/H] = -2.286 + 0.057 W a + 0.071 Wl (1) 

which was obtained by the least-square method. 

This calibration is valid for the regions of the 
(V — I, Mj, ECa) space covered by the calibra- 
tion clusters. It can be applied directly to the 
Fornax data below [Fe/H]~ —1.2 dex. However, 
Figure 5 shows that part of the Fornax data lies 
outside that region. There is no calibrating clus- 
ter that is simultaneously as bright, as blue and 
with as large ECa as the metal-richest Fornax tar- 
gets. The Mil sequence has similar colors and 
ECa but that cluster is not populous enough to 
contain stars as bright as in the Fornax sample. 
This means that the metallicity determination for 
the metal-richer part of our sample is not entirely 
constrained by the cluster calibrators and depends 
on a linear extrapolation of the ECa-Wo relation 
of about one magnitude in Mj. Until calibrators 
are available in that part of the parameter space, 
there will remain some uncertainties on the high- 
est [Fc/H] inferred from the Ca II triplet for the 
brightest red giants. 

The indications given about the extrapolation 
by the theoretical calculations of the Appendix 
are ambiguous. On the one hand, the theoreti- 
cal isometallicity sequences in the Mf-ECa plane 
show very little curvature between Mj = — 2 and 
Mj = —4. This would support a linear extrapola- 
tion to higher Mj. However, the theoretical rela- 
tions also indicate a steeper slope of the isometal- 
licity lines for higher mctallicities, whereas the ob- 
served cluster data is compatible with a constant 
slope. 

On the observational side, Ca II triplet data 
for red giants in the LMC (which are as bright 
and as blue as the Fornax stars with large Ca II 
equivalent widths) give a strong indication that 
the extrapolation of the calibration for high lu- 
minosities and high mctallicities is not linear in 
Mj. The LMC data provide complementary infor- 
mation to extend the calibration to stars younger 
and brighter than those of the calibrating clusters. 
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Recently, Cole et al. (2000) have collected Ca II 
data for an extensive sample of red giants in the 
LMC. The authors encountered the same difficul- 
ties in the high-metallicity end of their sample that 
we noted above in the case of Fornax. The prob- 
lem is even bigger for the LMC, because it is on 
average more metal-rich than Fornax. These au- 
thors dealt with the problem by taking M67 out of 
their calibration altogether, consequently extrapo- 
lating the W'-[Fe/H] relation in metallicity as well 
as the SCa-W relation in magnitude. The intro- 
duction of M67 in their calibration introduces a 
shift of the order of 0.4 dex in their most metal- 
rich objects (their Fig. 6). Fortunately there are 
other reliable constraints on the metallicity distri- 
bution of stars in the LMC, independent of the 
red giant Ca II measurements. Most studies indi- 
cate a typical value of about [Fe/H]~ —0.3 for the 
mean metallicity of the youngest populations (e.g., 

Luck et al. 1998 who obtain [Fe/H] 0.3 for the 

mean metallicity of Cepheids; the age-metallicity 
relation of Pagel & Tautvaisiene 1998 that cul- 
minates at [Fe/H]~ —0.2 based on observational 
data listed therein). The LMC Ca II data of Cole 
et al. (2000) are compared to the Fornax data 
on Figure 8. It shows that the LMC Ca II data 
does reach the luminosities of our Fornax sample, 
so that the better-constrained LMC data can be 
used to check the extension of the Ca II metallic- 
ity calibration to higher luminosities. 

We accordingly adapt our calibration for high 
luminosities by imposing a value of [Fe/H]= —0.3 
at Wo=5.5 A, corresponding to the metal-rich end 
of the LMC sample. The junction between the 
two calibrations is fixed at [Fe/H]=— 1.2 dex, the 
value where the Fornax red giants start having 
CMD locations different from the globular clus- 
ters. The resulting interpolation has very little 
curvature and is practically equivalent to the lin- 
ear formula: 

[Fe/H] = -2.87 + 0.47 W (2) 

This calibration is shown as a dotted line on 
Figure 7. Calibration (1) is valid up to Mj ~ 
—2, while the modified calibration (2) is applicable 
for Mi around —3 and brighter, like our Fornax 
sample. 

This procedure, which appears justified at this 
moment, should be checked by observing luminous 



red giants in open clusters that span wide ranges in 
age and composition, and more directly by observ- 
ing at high dispersion some of the metal-rich red 
giants in Fornax. We have examined the database 
of Cenarro et al. (2001b) for open cluster stars 
that overlap the Fornax stars in luminosity and 
colour. Only two stars in the cluster NGC 7789 
are as luminous as some of the Fornax stars, but 
these stars are significanlty redder. They do not, 
therefore, provide a direct test of the calibration 
technique taht we have applied to Fornax. Note 
that the discussion in the following sections does 
not depend strongly on the exact values of the 
highest metallicities in Fornax, the crucial factor 
being the presence in the sample of a large number 
of stars with very high values of SCa. 

To apply the Mj calibration to the Fornax data 
we use E(V— 1)=0.07 for Fornax, which is the mean 
value that Buonanno et al. (1998) derived for 
four globular clusters in Fornax 9 , V(HB)=21.28, 
the mean V magnitude of the horizontal branch of 
the same 4 clusters, and My (HB) =0.47, from the 
models of Demarque et al. (2000). These numbers 
give (m-M) v =20.81 and (m-M)/=20.74. 

The resulting values of [Fe/H] for our sample 
are given in Table 2, and the resulting metallic- 
ity distribution is displayed in Figure 9. Figure 6 
plots the [Fe/H] as a function of I magnitude. As 
a reminder of the fact that the uncertainty on the 
calibration is higher for the metal-rich part of the 
Fornax sample, the values above [Fe/H]=— 0.9 are 
rounded off and indicated by a colon (":") in Ta- 
ble 2, and the corresponding bins are shaded in 
Fig 9. 

This metallicity distribution is broadly consis- 
tent with that of Tolstoy et al. (2001): in both 
cases, most of the stars in the sample have values 
in the range — 1.5 < [Fe/H] < — 0.7 and, if we con- 
sider Poisson error bars in the number of stars with 
metallicities above and below [Fe/H] = —1, the 
relative number of stars may be compatible. Our 
sample is, however, centered at significantly higher 
abundances (67% of stars with [Fe/H] > — 1, com- 
pared to 53% in the Tolstoy et al. sample, and the 
peak of the distribution is [Fe/H] = —0.9 com- 

9 The Schlegel, Finkbeiner & Davis (1998) reddening maps 
show that the variations of E(V— I) accross the field of For- 
nax are at most a few hundredths of magnitude, negligible 
for our purpose. 
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pared to [Fe/H] = -1.2 in Tolstoy et al.) These 
differences may be in part explained by the pop- 
ulation selection effects resulting from the differ- 
ent brightness of the two samples: our sample is 
within ~1 mag below the tip of the RGB, while 
the Tolstoy et al. sample is between 1 and 3 mags 
fainter than the tip. As we will show in the next 
section, selecting bright red giants may imply a 
bias towards young ages and, considering a gen- 
eral metallicity law in which metallicity increases 
with time, higher metallicities. Our results are, 
in any case, statistically more solid, since we have 
more than three times the number of stars of Tol- 
stoy et al. (2001). Finally, note the much higher 
quality of our spectroscopic data, reflected in our 
much lower error bars that those inferred by Tol- 
stoy et al. (0.23 dex on average vs. 0.09 dex of 
mean error in our measurements). The difference 
in the quality of the spectra is also evident if one 
compares their hgure 10 with our figure 4. In their 
figure 10, Tolstoy et al. show a good spectrum 
of Fornax, with S/N~30, which is comparable in 
quality to all our spectra. The star to which this 
spectrum belongs, however, is about 2 magnitudes 
brighter than the rest of their Fornax sample. 

3.2. Inferred metallicity distribution for 
the whole population 

The metallicity distribution of our red giant 
sample is an indirect reflection of the metallicity 
distribution of Fornax as a whole, since the bright- 
est red giants are only a specific sub-set of the 
galaxy. To infer the likely metallicity distribution 
of the underlying population in Fornax from our 
data, we computed synthetic CMDs and examined 
the relation between the "observed" and "real" 
metallicity distribution in these models. The mod- 
els use the synthetic CMD code ZVAR (described 
in Gallart et al. 1999) with the Bertelli et al. 
(1994) generation of Padova stellar evolutionary 
models, an IMF from Kroupa, Tout & Gilmore 
(1993) and an input star formation rate and AMR 
to produce a synthetic population. For a descrip- 
tion of these parameters and the way they are 
used, the reader is referred to Gallart et al. (1999). 
We used a constant star formation rate and exper- 
imented with several AMR. All AMR compatible 
with our Fornax data produce essentially similar 
results, and we give below the results for an AMR 
of z=0.0004 at agc=15 Gyr, z=0.002 at agc=1.5 



Gyr, z=0.008 at age=0 Gyr, with linear interpola- 
tion in between (the maximum age for the Padova 
isochrones that we use is 15 Gyr). On this syn- 
thetic population, we then apply our selection cri- 
teria (see Section 2.1) and compute how many ob- 
jects fall in different age and metallicity bins. The 
results are given in Tables 3 and 4. As the input 
star formation rate is constant, the distribution of 
the ages of the selected stars indicates the age bias 
relative to the underlying population (or, more 
precisely, to the population of stars with tu f e > r, 
such as K-dwarfs) i.e. the synthetic CMDs in- 
dicate that we selected preferentially stars with 
ages < 4 Gyr, by about a factor of a 30%. Con- 
sequently, the metallicity distribution of the RGB 
sample is also different from that of the underly- 
ing population. The values in Tables 3 and 4 are 
used in Section 4 to infer the age and metallic- 
ity distributions of the whole population from the 
distributions of our sample. 

3.3. Potential sources of uncertainty 

i) Possible differences in [Ca/Fe] between 
Fornax and the calibrating clusters. The 

equivalent width of the Ca II triplet depends on 
the abundance of Calcium, an a-element. Relat- 
ing it to the [Fc/H] scale therefore requires knowl- 
edge of the run of [Ca/Fe] with [Fe/H]. Among 
the calibrating clusters, there is significant vari- 
ation in [Ca/Fe] that is not uniquely correlated 
with [Fe/H]. The globular clusters that make up 
the calibration below the metallicity of 47 Tuc 
QFe/H]=-0.7) have on average [Ca/Fe]~ +0.3 
(Carney 1996; Sneden et al. 1997; Gonzalez & 
Wallerstein 1998; Shetronc & Keane 2000). Three 
studies of 47 Tuc itself (see Carney's review) in- 
dicate that it has [Ca/Fe] ~ 0.0, as does the open 
cluster M67 (Tautvaisiene et al. 2000). The two 
globular clusters that are more metal-rich than 47 
Tuc, NGC 6528 and NGC 6553 have [Ca/Fe] ~ 
+0.25, according to Carretta et al. (2001). This 
result and the measurements of [Fe/H] by these au- 
thors yield [Ca/H] values that actually exceed the 
value obtained by Tautvaisiene et al. (2000) for 
M67, which is surprising because these clusters lie 
below the M67 sequence in the M/-SCa diagram 
(see Fig. 5). While this apparent inconsistency 
may be simply due to observational error, it is im- 
portant to note that the M/-ECa diagram is not 
only sensitive to [Ca/H] but also to the values of 
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Te and log g of the stars evolving on the RGB (see 
Appendix) , which are in turn complex functions of 
the overall metallicities of the stars and their ages. 
It is not certain that the values of Wo that are de- 
rived from the M/-£Ca diagram should be more 
closely correlated with [Ca/H] than with [Fe/H]. 
Further theoretical modeling of the Mj-£Ca dia- 
gram is urgently needed to clarify this point and to 
put the question of its age sensitivity on a firmer 
footing. 

The possibility remains that the Fornax stars 
have so different [Ca/Fe] than our calibrating ob- 
jects that their positions in the Mj-£Ca diagram 
in Fig. 5 are not indicative of their metallicities. 
However, the very recent high dispersion analy- 
ses of three red giants in Fornax by Shetrone et 
al. (2003) suggest that this is not the case. They 
obtained [Ca/Fe] = +0.23, +0.21, and +0.23 for 
stars that have [Fe/H]= -1.60, -1.21 and -0.67, 
respectively. Since these values of [Ca/Fe] are near 
the mean of our calibrating clusters, we highly 
doubt that the Fornax stars in our sample have 
so different [Ca/Fe] values that the Mj-£Ca dia- 
gram produces for them an erroneous ranking by 
[Fe/H] . 

ii) AGB stars. Synthetic population simula- 
tions (see Section 3.2) indicate that about a quar- 
ter of the stars in our sample are probably AGB 
stars rather than RGB stars, and that they are 
brighter on average. However, both Cole et al 
(2000) and our simulation described in the Ap- 
pendix lead to the conclusion that metallicities de- 
rived for AGB stars with the Call triplet method 
suffer from very small systematic biases compared 
to RGB stars. This is also confirmed by the fact 
that there is no systematic trend of metallicity 
with magnitude in our sample, as would be the 
case if AGB stars -brighter on average- were pro- 
ducing biased metallicities. 

iii) Possible field stars contamination. It 
is very unlikely that the high-SCa stars in our 
sample are field stars, which are expected to be 
primarily dwarf stars of the thin and thick disk 
populations. From the Ca II measurements of 
Cenarro et al. (2001a), we estimate that even 
the most metal-rich of these high-gravity stars will 
have SCa< 5.4 A, and therefore cannot be con- 
fused with the Fornax stars with SCa > 6 A. The 
red giants of these populations are too bright to 
be confused with our sample of Fornax stars, and 



the small number of halo red giants along the line 
of sight to Fornax is unlikely to contain any with 
SCa larger than the 47 Tuc values. 

There may be a few field objects in the low-SCa 
portion of the Fornax sample. Indeed some of the 
blue, low-SCa objects in our sample are situated 
on the blue side of the Fornax RGB, where the 
density of the field objects is similar to that of 
the Fornax objects. Their color, magnitude and 
SCa indicate that they can be intermediate-age, 
low-metallicity Fornax red giants but also metal- 
rich thick disc dwarfs. However, the radial velocity 
data indicates that any such field contamination is 
small. 

iv) Uncertainty in the distance modulus. 

The uncertainty in the distance modulus affects 
the calculation of Mj in the metallicity calibra- 
tion. With our calibration, a change in ±0.12 mag 
on the distance modulus of Fornax, corresponding 
to the uncertainty of Saviane et al. (2000), modi- 
fies the recovered [Fe/H] by ^0.03 dex. This is a 
negligible value in view of the other uncertainties. 

4. Analysis: the star formation and chem- 
ical enrichment history of Fornax 

Although individual values of [Fe/H] at the 
higher end have to be taken with caution, the Ca II 
triplet data put very strong lower limits on the 
possible metallicity and upper limit on the pos- 
sible age for the majority of the Fornax sample. 
This is a very important result of our study, that 
in turn strongly constrains the possible chemical 
enrichment and star formation histories for For- 
nax. In this section, we compare the metallicities 
of the stars with their positions in the CMD, and 
further examine the age-metallicity relation that 
we have derived in Fornax. 

4.1. Compatible chemical enrichment his- 
tory 

It is very instructive to compare the computed 
metallicities with the positions of the stars in the 
CMD. In the upper left panel of Fig. 5, the Fornax 
stars are plotted along with the RGB's of some of 
the clusters for which Ca II triplet data are avail- 
able. The majority of the Fornax stars lie between 
the RGB's of M15 and NGC 1851. If all of these 
stars were very old, comparable in age to the glob- 
ular clusters, then their metallicities should lie be- 
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tween [Fe/H]=— 2 and —1.1, and their Ca II equiv- 
alent widths should lie between the sequences for 
these clusters. The positions of the red giants in 
the metal-rich and young open cluster Mil are 
also plotted in this diagram to underscore the fact 
that a metal-rich but young population is indis- 
tinguishable on the RGB from a metal-poor, old 
population. 

In the Mi - SCa plane (see upper right panel 
of Fig. 5), most of the Fornax stars lie above the 
line for NGC 1851, which indicates that a large 
fraction of the Fornax sample is more metal-rich 
than this cluster. Because these same stars are 
bluer than the RGB of NGC 1851 in the CMD 
(Fig. 5), they must be substantially younger stars. 
This age effect is very evident in the case of the 
Mil stars, whose metallicities and ages are well 
documented. 

Figure 10 illustrates the fact that the large SCa 
values for the Fornax stars combined with their 
blue color in the CMD impose strict upper limits 
on their ages. Figure 10 indicates the relation be- 
tween the predicted color of the RGB of a cluster 
and its age, for different metallicities, according 
to the Yonsei-Yale (Yi et al. 2001) isochrones. 
The red limit of the Fornax RGB is indicated as 
a dashed line. For stars that are of metallicities 
equal or higher than [Fe/H]=— 0.68 for instance - 
which is the case of many of the Fornax sample 
stars according to the Ca II triplet data- ages 
have to be lower than 2 Gyr in order for the red 
giant sequence to be as blue as observed. There- 
fore, our Ca II triplet data implies not only un- 
expectedly high abundances for the Fornax stars, 
but also quite young ages. 

Two immediate, qualitative conclusions from 
our data are the following: First, the AMR of 
Fornax must be rather tight, otherwise the scatter 
would produce older metal-rich stars, that would 
be visible near to the 47 Tuc locus and are not 
observed. Second, the stars of metallicities near 
[Fe/H]= —1 must have ages below 5 Gyr approxi- 
mately, while those near [Fe/H]= —0.6 must have 
ages below ~2 Gyr, in order to be at the position 
observed in the CMD. 

Age estimates were computed for our objects 
by calibrating the age dependence of the RGB 
with Padova stellar evolution models (Girardi et 
al. 2000). Figure 11 shows the resulting age es- 
timates versus the metallicities obtained from the 



Ca II triplet. The horizontal bars show the 95% 
confidence interval of the age determination given 
the observational uncertainties. The most likely 
age values and the confidence intervals were com- 
puted following the Bayesian approach of Pont & 
Eyer (in prep.). Bayesian probability theory al- 
lows a realistic computation of the probability dis- 
tribution function of a derived quantity, here the 
age, when the relative uncertainties are high (see 
for instance Sivia 1996). In such cases, a sim- 
pler approach often leads to systematic biases and 
underestimated error intervals. The possible pres- 
ence of AGB stars would slightly increase the com- 
puted ages for a fraction of the sample (by ~0-l 
Gyr for young objects), the AGB being slightly 
bluer than the RGB in most cases. Note that 
reasonable changes in the mctallicity calibration 
do not affect the fundamental features of these re- 
sults, young ages being imposed to a large fraction 
of the sample by the fact that the metal-rich For- 
nax giants are much bluer than globular clusters 
of similar metallicities. 

Fig. 11 gives an indication of the age-metallicity 
relation of Fornax. The constraint on ages get 
stronger for stars with ages below 5 Gyr. The 
emerging picture (ensuring coherence between the 
spectroscopic and photometric data on Fornax red 
giants) is that Fornax underwent an initial stage 
of enrichment reaching [Fe/H]~ —1 about 3 Gyr 
ago. Then, sustained star formation progressively 
increased the ISM [Fe/H] to a recent value of at 
least ~ —0.5 dcx. The metal enrichment curve 
of a closed-box model is plotted on Fig. 11 for 
comparison. 

Table 5 shows the result of sorting the ages of 
the stars into four intervals and provides an indi- 
cation of their age distribution. To infer informa- 
tion about the star formation rate of Fornax as 
a whole from this distribution, the statistical cor- 
rections calculated in Section 3.2 from synthetic 
models were applied. As a result, Table 5 would 
imply a significant increase in the star formation 
rate in the last ~4 Gyr. 

These conclusions are affected by the uncertain- 
ties in the relative lifetimes of stars of different 
ages near the tip of the RGB. In fact, our ex- 
perience with synthetic CMD studies (Gallart et 
al. 1999; Gallart, Aparicio & Bertelli 2003) shows 
that the number counts of stars near the tip of 
the RGB is difficult to match consistently with the 
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number of stars in other parts of the CMD. Zoccali 
& Piotto (1999) found good agreement between 
the observed and theoretical luminosity functions 
of the main-sequence, subgiant-branch and RGB 
of old globular clusters in a range of metallicities. 
A similar study involving clusters of different ages 
is highly needed. 

4.2. Other clues from the CMD 

Is our picture compatible with parts of the 
CMD other than the RGB? The Fornax CMD con- 
tains numerous features associated with a partic- 
ular age and metallicity, namely a weak horizontal 
branch, a prominent red clump with a long tail at 
the bright end, and a main sequence extending to 
Mj ~ —0.5. These features indicate respectively a 
weak old metal-poor population, an intermediate- 
age population of various ages and metallicities, 
and a very recent population, down to ~500 Myr 
old. All of these are qualitatively compatible with 
the picture obtained from the RGB and the cal- 
cium triplet. 

To check the compatibility of our AMR with the 
observed CMD more quantitatively, we produced 
a synthetic CMD using the AMR obtained above 
-approximated by a linear interpolation between 
the values (Gyr, Z)= (15, 0.0002), (13, 0.0008), 
(2, 0.0036), (0,0.008)- and constant star formation 
rate between 12 and 0.5 Gyr (see Figure 12, and 
its caption, for more details). In Figure 12, this 
synthetic CMD (middle panel) is compared to the 
Fornax observations (left panel). A visual compar- 
ison shows that the main features are well repro- 
duced in the model, except for the fact that the po- 
sition of the RGB is not perfectly matched, which 
is a well-known problem of current stellar evolu- 
tionary models Since the main sequence is the the- 
oretically best understood feature in a CMD, par- 
ticularly important is the good agreement between 
the position and shape of the observed and syn- 
thetic upper main-sequence, which are both well 
centered within the corresponding reference lines. 
Since the AMR has been obtained from informa- 
tion on the RGB exclusively, this is a completely 
independent consistency check indicating that the 
spectroscopically derived metallicities are in good 
agreement with the observed CMD. 

As a comparison, we computed a synthetic 
CMD with a slightly different AMR (as can be 
deduced from Saviane et al. 2000) whose main 



difference with the AMR inferred by us is a lower 
metallicity for the 1-2 Gyr old population. This 
CMD is represented in the right panel of Fig- 
ure 12. Note that the post-main sequence features 
are quite similar in the two synthetic models, and 
that the main difference is in the position of the 
upper main sequence, which is too blue when the 
lower metallicity AMR is used. This further indi- 
cates that the high metallicity, young end of the 
AMR is a feature not only suggested by the spec- 
troscopic measurements, but necessary to obtain 
a good match with the young part of the observed 
CMD. 

5. Discussion: Fornax as a Local Group 
dwarf galaxy 

The metallicity distribution and AMR that we 
derive for Fornax make it a complex system, more 
reminiscent of the LMC (Pagel & Tautvaisiene 
1998) or the Milky Way's disk (Rocha-Pinto et 
al. 1996, Edvardsson et al. 1993) than of other 
dwarf spheroidals. Like those two systems and 
unlike most other dwarf spheroidals, Fornax has a 
rather metal-rich main population, with a tail of 
old and metal-poor stars. Also like these two sys- 
tems, its AMR is compatible with a rapid initial 
enrichment, followed by a period of slower enrich- 
ment, and then a recent acceleration of the enrich- 
ment. 

Fig. 13 compares the cumulative metallicity dis- 
tribution of Fornax with that of a "closed-box" 
model and the Galaxy. Although there is also a 
departure from the closed-box model in Fornax, 
the distribution is significantly nearer to a closed- 
box relation than the Galaxy. This is probably 
even more true when the outer parts of Fornax 
are taken into account 10 . The inclusion of the 
outer parts in the cumulative metallicity distri- 

'Fornax is known to be spatially segregated, with a larger 
fraction of intermediate-age and young stars in the inner 
parts (Saviane et al. 2000). Our results, therefore, are rep- 
resentative only of the area of the galaxy that they cover, 
the central 11' minutes in radius of Fornax. This radius is 
near the exponential radius given by Mateo (1998) for For- 
nax. Varying assumptions about the density profile of For- 
nax within reasonable limits, the portion of the galaxy sam- 
pled by our data ranges between 14% (flat profile out to 32 
arcmin, then exponential profile with 7 = — 0.1 arcmin - 1 
(Demers et al. 1994) out to total radius 71 arcmin (Ma- 
teo 1998)) and 24% (exponential profile throughout with 
7 = — 0.1 arcmin -1 ). 
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bution will proportionally increase the number of 
very metal-poor objects. 

Consequently, there appears to be no need to 
invoke a large infall of unenriched gas on the 
galaxy or a large loss of enriched gas. Moreover, 
the almost complete absence of young, metal-poor 
stars would seem to exclude a scenario of accre- 
tion of relatively pristine gas for late star forma- 
tion episodes. 

5.1. Where is the gas? 

The Fornax results strongly reinforce the mys- 
tery posed by the complete absence of gas in 
dwarf spheroidal galaxies (Knapp, Kerr & Bowers 
1978). Fornax seems to have been producing large 
amounts of stars continuously until very recent 
times. Then, as we observe it, it has stopped form- 
ing stars and disposed of every trace of gas. We 
could be observing Fornax just after its last star 
formation episode, similar to Carina and Leo I. 
Dwarf irregulars could be equivalent galaxies form- 
ing stars at the present time. However, the situa- 
tion would still be statistically unlikely. If Fornax 
has been forming stars up to 200 Myr ago (Saviane 
ct al. 2000), it would be a surprising coincidence 
to observe it this close to the time when its gas was 
rapidly enriched and then completely expelled. 

A more satisfactory solution would be if some 
gas is hidden but still at the disposal of For- 
nax. Blitz & Robishaw (2000) reported several 
dSph galaxies with HI in the vincinity. While sev- 
eral of these detections were not confirmed (e.g., 
Young 1999), and some even retracted (Blitz, priv. 
comm.), the HI around the Sculptor galaxy has 
been confirmed by Carignan ct al. 1998. Fornax 
was unfortunately not included in the Blitz & Ro- 
bishaw study because it was outside the area cov- 
ered by the Dwingaloo survey. Given its history 
of continuous and recent star formation, we would 
predict that significant amounts of gas associated 
with Fornax might still be detected in its vicinity. 

Another conclusion of Blitz & Robishaw (2000) 
is that dSph nearer than 250 pc to a giant galaxy 
(the Galaxy or M31) contain much less gas, an 
effect that the authors attribute to ram-pressure 
stripping. Fornax is well within this distance lim- 
its, and we would accordingly expect it to contain 
little or no gas. We may still, then, be viewing For- 
nax by accident just after its "gas death", during 



what may be its first passage very near the Milky 
Way. The precise determination of the proper mo- 
tion of Fornax and a computation of its likely orbit 
will provide a test of this scenario (see Piatek et 
al. 2002 for a preliminary estimate). 

5.2. Comparison with evolutionary mod- 
els of dwarf galaxies. 

It is interesting to compare the properties of 
Fornax with the predictions by Mac Low & Fer- 
rara (1999, MF99) and Ferrara & Tolstoy (2000, 
FT00) . In order to identify where is Fornax in the 
MF99 parameter space, we will derive its baryonic 
mass, and an estimate of its likely SNe rate from 
data in the literature. The total (dark plus bary- 
onic) mass of Fornax is — 7 x 10 7 (Mateo 1998); by 
assuming the (M/L) to t,v = 4.4 given in the same 
paper, and a baryonic M/L ~ 1.0 (calculated for a 
population with constant star formation rate from 
12 to 0.5 Gyr ago), this corresponds to a baryonic 
mass of M g = 1.6 x 10 7 M Q . For a constant star 
formation rate, it corresponds to a rate of 1400 
M Q Myr" 1 . For a Salpeter IMF, and this star 
formation rate, 170 M Q Myr" 1 of stars would be 
produced in the mass interval of SNe II progeni- 
tors, 10-100 M Q , with mean mass ~ 20M Q . These 
figures imply a rate of ~ 8.5 SNe/Myr, or a SNe 
luminosity ~ 2.7 x 10 38 ergs s _1 . If we input these 
data in Tables 2 and 3 of MF99, we find that the 
mass ejection efficiency expected for Fornax would 
be modest, ~ 1 x 10" 2 , while all the metals should 
be ejected. The fact that there is a substantial 
metal enrichment in Fornax seems to be in con- 
tradiction with these models, and would indicate 
that the metal ejection efficiency must be in some 
way overestimated in them. 

Is the model more representative for the small- 
est dSphs? Making the same calculation as above 
for a system with baryonic mass M g = 3 x 10 5 M Q , 
and assuming a constant star formation rate ac- 
tive in them for a maximum of ~ bGyr, we get an 
upper limit for the star formation rate of 60 M 
Myr" 1 , or 0.3 SNe Myr" 1 . This is the lower limit 
of SNe rate considered by Mac Low & Ferrara, 
who predict that these galaxies are also expected 
to keep some gas through their lifetimes but to lose 
all metals produced by type II SNe. This may be 
compatible with their enrichment being substan- 
tially lower than that of Fornax. However, lower 
star formation rates and therefore lower SNe rates 



14 



than calculated above must be invoked for some 
fraction of the galaxy's lifetime in order to explain 
the still noticeable enrichment of the gas reflected 
in the metallicity of their stars (Shetronc, Cote & 
Sargent 2001 and several earlier studies). Tamura 
ct al. (2001) have proposed a possible solution 
that may work at the lowest metallicities. They 
note that star formation may be strongly regulated 
in low metallicity gas by the lack of coolants (Nishi 
& Tashiro 2000) . But once the metallicity rises to 
[Fe/H]~ —2, the metals become effective coolants, 
and the subsequent star formation should produce 
SNe that expel their metals as described by MF99 
and FT00. This is inconsistent with the observa- 
tional data for Draco, Sextans, and Ursa Minor 
(Shctrone et al. 2001), which show that they con- 
tains substantial numbers of stars that are more 
metal-rich than this limit. 

One possible reason for the apparent overes- 
timate of metal ejection efficiency in the MF99 
and FT00 models is the fact that they consider 
that all the mass injection occurs in the central 
100 pc of the galaxy. As MF99 discuss, several 
small bubbles distributed in a larger area must 
be less effective in blowing away the ISM than 
a single central cluster. The core radius of For- 
nax is 460 pc and that of the smaller galaxies is 
~ 200 pc. Since dwarf irregular galaxies show spa- 
tially and chronologically segregated sites of star 
formation (as beautifully shown by Dohm-Palmcr 
ct al. 2002), it is possible that the assumption of 
one large central site for the dSph galaxies is also 
incorrect. 

Further hints may be provided by high reso- 
lution spectra abundances of individual elements. 
Shetrone et al. (2001, 2003) note that all the dSph 
have an under-abundance of even-Z elements com- 
pared with the Milky Way halo, which may indi- 
cate a large contribution by type la SNe, cither be- 
cause SNe II are less abundant in these galaxies or 
because their ejecta are not retained. Could it be 
that the ejecta from type la SNe are retained more 
than that of type II, and that this provide some 
extra enrichment? One property of SNe la which 
is different from SNe II is that they are expected 
to be isolated events, i.e. less spatially clustered 
and likely to be less synchronized in time than 
SNe II. According to FT00, isolated SNe events 
would not cause blowout and would have a lower 
metal ejection efficiency. There seems to be one 



difference between the Fornax abundance patterns 
and that of other dSph, according to Shetrone et 
al. (2003): while there is a general trend of de- 
creasing [a/Fc] abundance with increasing metal- 
licity in the observed dSph, in Fornax the even-Z 
ratios appear flat to slightly rising, even if still 
under- abundant. All this seems to indicate that 
the enrichment mechanisms are similar in dSph 
regardless of their mass. The slightly different 
tendency of even-Z ratios with metallicity in For- 
nax may simply reflect the more extended period 
of star formation in this galaxy compared to the 
other dSph. This discussion, however, is based on 
very few stars. High resolution spectra of a large 
number of stars will provide very important infor- 
mation to understand the enrichment mechanisms 
in dSph galaxies, and the processes driving their 
evolution. 

6. Summary and conclusions 

We obtained spectroscopy in the Ca II triplet 
region of the spectrum for 117 stars in the Fornax 
dSph galaxy using FORS1 at the VLT. We de- 
rived metallicities for them using our own SCa-M/- 
[Fe/H] calibrations, based on observations of Ca II 
line strengths in 11 globular clusters plus M67, 
Mil and the LMC. We show evidence that a cali- 
bration against M/ instead of the most commonly 
used one in terms of V-V(HB) may produce more 
accurate results for the young populations present 
in Fornax. Our main conclusion is that there is a 
tail of very metal-rich RGB stars in Fornax, with 
metallicities ranging from approximately —0.7 to 
—0.4, whose existence could not have been sus- 
pected from their color in the CMD. The metallic- 
ities of these stars, combined with their relatively 
blue colors, imply that they have ages of the order 
of 2 Gyr. 

We use synthetic CMDs to examine how the 
age and metallicity distribution of our sample is 
related to that of the whole population. Combin- 
ing the metallicity data with constraints from the 
position of the same stars in the CMD, we find that 
the Ca II data are compatible with only a narrow 
range of possible AMR, which in turn, allow to put 
some constraints on the SFH. It appears that For- 
nax underwent very substantial chemical enrich- 
ment in the last few Gyr, from [Fe/H]~ —1.0 up 
to [Fe/H]~ -0.4. Its AMR is nearer to a Closed- 
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Box model than that of the Galaxy. The large 
number of young (< 4 Gyr old) stars that we infer 
from the combination of their metallicity and their 
position on the CMD would imply a star forma- 
tion rate somewhat rising (by a factor of ~2) in 
the last 2-4 Gyr. However, this last result depends 
on an extrapolation from the RGB population to 
the total population that is not very reliable with 
the present models. 

This scenario corresponds remarkably well with 
observed features in the CMD other than the 
RGB, and in particular, with a robust feature (in 
terms of our stellar evolution knowledge) as is the 
main sequence. Indeed, the derived AMR, and 
in particular, the high metallicity tail, produces 
synthetic CMDs with upper main-sequence colors 
totally compatible with the observed one, whereas 
a lower metallicity for the youngest stars would 
produce too blue a main sequence. 
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A. Theoretical sensitivity of the Ca II triplet to age and metallicity. 

Under the usual assumption that abundance, effective temperature (T e ) and surface gravity (log g) are 
the main parameters determining the strengths of absorption lines, the combined equivalent width of the 
Ca II triplet (SCa), evolves along a surface in the (T e , log g, [Ca/H], SCa) space. Consequently, to relate 
abundances to the measurable parameter, ECa, the parameters T e and log g must be fixed in some way. It 
was found empirically (Da Costa 1998; Armandroff & Da Costa 1991; Olszewski et al. 1991) that for globular 
clusters a plot of SCa against absolute V or I magnitude or against the closely related quantity V-V(HB), the 
difference in V magnitude between a red giant and the mean V of the horizontal branch, produces an accurate 
ranking of globular clusters by metallicity. The precision in [Fe/H] that one can obtain for individual red 
giants using ECa and, for example, V-V(HB) is quite high (< ±0.2 dex), and this technique has developed 
into the most popular way of using low resolution spectra to estimate the abundances of stars in globular 
clusters and dSph galaxies. Before applying this method to the red giants in Fornax, we need to understand 
why it works well for globular cluster stars and its limitations for a mixture of ages and compositions. 

The globular clusters in the Milky Way are very old, and for the purposes of our discussion their age 
differences can be ignored. With the exception of a very few clusters, most notably ui Cen, the stars 
within a globular cluster have very nearly the same chemical composition. These facts explain why in the 
CMD, the RGBs of globular clusters are tight sequences of stars that are offset from one another in order 
of the metallicities of the clusters, as predicted by stellar evolutionary theory for clusters of uniform age 
but varying composition. Since T e , log g, and luminosity vary monotonically along the RGB, one can use 
measures of T e , typically B-V or V-I, or measures of luminosity, My, Mj, or V-V(HB), which is independent 
of both reddening and distance, as indicators of position on the RGB. It is known empirically that plots 
of SCa against any one of these luminosity indicators provides better metallicity discrimination than a plot 
against an indicator of T e . To understand why this is so let us consider the sensitivity of SCa to T e and 
log g using the calculations that Jorgensen, Carlsson, & Johnson (1992) made of SCa for stars in the ranges 
4000 <T e < 6600, 4.0 < logs- < 0.0, and -1.0 < [A/H] < +0.2, under the assumption of NLTE. For 
estimates of T e and \ogg for the red giants, we will use the Padova isochrones (Girardi et al. 2000). It is 
important to note that the calculations by Jorgensen et al. (1992) and by Girardi et al. (2000) assumed 
solar mixes of elements. While there is clear observational evidence that [Ca/Fe] is not solar in at least some 
Fornax stars (Shetrone et al. 2003), their departures are small in comparison to the effects that we describe 
below, which provide only a qualitative description of the behavior of the Ca II lines. To avoid confusion with 
predictions for specific abundance ratios, we will use [A/H] to indicate overall metallicity when discussing 
the Call line strengths. 

Figure 14 shows the most luminous portions of the RGB in the log T e - log g plane for four different 
metallicities (Z=0.001, 0.004, 0.008, and 0.019, which assuming Z© = 0.02 corresponds to [Fe/H]= -1.30, 
—0.70, —0.40, and —0.02, respectively). The quantities are plotted so that this diagram resembles an H-R 
diagram, in that T e increases to the left and luminosity increases toward the top. Representative lines of 
constant I and V absolute magnitude are shown. As in the H-R diagram, the RGB's of different Z form 
a nearly parallel sequence of curves. At every log g, the most metal-rich RGB has the coolest T e . A 
vertical line in Figure 14 indicates that there is a T e below which the Call lines are seriously blended with 
a TiO band. The location of this line was estimated from the V-I colors of red giants having significant 
TiO absorption in the very old, metal-rich open cluster NGC 6791 (Garnavich et al. 1994) and the more 
metal-poor globular cluster 47 Tuc (Lloyd Evans 1974). While this line is approximate, it illustrates the 
observed fact that ECacannot be measured to the RGB tip in 47 Tuc ([Fe/H] = —0.7) and in more metal-rich 
globular clusters. 

According to the calculations of Jorgensen et al. (1992), at fixed T e and log g, SCa increases with 
increasing [A/H]. At fixed [A/H] and log g and over the T e range of the RGB's in Figure 14, SCa increases 
with increasing T e , with the exception that for [A/H]=— 1, it decreases slightly for log g > 2. At fixed 
[A/H] and T e , these calculations show that for log g < 2, SCa increases with decreasing log g. This is much 
more pronounced at [A/H]=0.0 than at lower abundances, and in fact for [A/H] = — 1 SCa has little gravity 
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sensitivity except between log g = 1 and 0. Figure 14 illustrates that going along a line of constant absolute 
I or V magnitude from the Z=0.001 RGB to the one for Z=0.019, both log T e and log g decrease. The 
decrease in T e tends to weaken ECa, while the decrease in log g tends to strengthen it. The two effects 
roughly cancel. The dependence of ECa on [A/H] ensures that RGB's of different Z will be separated in a 
plot of ECa against absolute I or V magnitude or V-V(HB), as is observed. 

A similar argument explains why ECa changes by only a small amount along the RGB of a globular 
cluster. As luminosity increases, the accompanying decreases in T e and log g weaken and strengthen ECa, 
respectively, which results in only a modest increase in ECa with luminosity (5T,Ca,/5Mj ~ 0.5). This explains 
why in the absolute magnitude vs. ECa diagram, the RGB's of globular clusters are well approximated by 
straight lines. It is remarkable, however, that these lines have the same slope to within the errors over a 
wide range in mctallicity. In the log T e - log g plane, the RGB's of different Z are almost straight lines and 
approximately parallel (see Fig. 14). If the partial derivatives of SCa with respect to log T e and log g are 
nearly constant for all mctallicitics, then the RGB's will be lines of similar slope in the absolute magnitude 
vs. ECa diagram. This is not the case according to the calculations of Jorgensen et al. (1992), which predict 
that Ca II strength increases more rapidly with decreasing log g at high metallicity than at low, which should 
cause an increase in the slope with increasing metallicity in the absolute magnitude vs. ECa diagram (see 
Fig. 16). It is surprising that this effect is not clearly evident in the data (see Section 3.1). 

We have also illustrated in Figure 14 the curve for a constant color of V— 1=1.4, which as expected is 
nearly a line of constant T e . This curve illustrates why plotting ECa against V— I or another color sensitive 
to T e produces relatively poor separations between the RGB's of different metallicities. Along a line of 
constant V— I, T e remains essentially constant, but log g increases substantially from the metal-poor to the 
metal-rich RGB. These differences in gravity weaken the dependence of ECa on metallicity and tend to push 
together the RGB's of globular clusters in plots of W(Ca II) against color. This effect and the sensitivity 
of color to errors in the reddening corrections make plots of ECa against color much less useful than plots 
against absolute magnitude. 

In this article we calibrate ECa against metallicity for a mixed-age population. It is very instructive 
to make similar comparisons of RGB's holding Z fixed, but varying the ages of the populations. We have 
plotted in Figure 15 the Padova isochrones for 12.6 and 0.25 Gyrs and Z=0.019. A line of constant absolute 
I magnitude is also plotted, and along this line we have marked the locations of the isochrones for this 
composition and ages of 12.6, 4.0, 1.0, 0.5, and 0.25 Gyrs. For ages < 1 Gyr, the RGB terminates near 
Mj —3, and the more luminous stars of these ages are AGB stars. Since the AGB and RGB are almost 
coincident in the log T e - log g plane (see Fig. 15), the ambiguity about whether a star belongs to the AGB 
or the RGB has only a small effect on the interpretation of its ECa. Figure 15 illustrates the important point 
that at a constant absolute magnitude, the isochrones for younger ages have higher T e and larger logg than 
do ones for older ages. Because ECa is affected in opposite ways by these differences, age differences have 
only small effects on plots of ECa against Mj or My. The differences in T e and log g between the 12.6 
and 4.0 isochrones are so small that their effects on ECa nearly perfectly cancel. This is not the case for 
differences between the 12.6 isochrone and the ones for ages < 1 Gyr, but nonetheless the Mj-ECa remains 
a very useful diagnostic of metallicity even if there is a large range in age in a stellar population. 

Fig. 16 illustrates the weak sensitivity of the Mj-ECa to age differences. We have used the approximation 
formulae that J0rgensen et al. (1992, eqs. 7, 9, & 11) derived from fits to their calculations of the sum 
of the equivalent widths of the two strongest Ca II lines (8542 and 8662A) and the Padova isochrones to 
compute ECa for stellar populations of different ages and compositions. The values of T e and log g that are 
required as input to the Jorgensen et al. (1992) relations for compositions [A/H]=0.0, —0.5, and —1.0, were 
obtained from the Padova isochrones for Z=0.019, 0.008, and 0.004, respectively. The last two compositions 
are not exact matches, but repeating the calculations using instead Z=0.004 and 0.001 for [A/H]=— 0.5 and 
— 1.0, respectively, showed that this has a negligible effect on the ECa vs. M/ diagram. It has only a small 
systematic effect on the ECa vs. V-I diagram in Figure 16 in the sense that the lines of constant composition 
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are separated slightly more in SCa than they would be if the correct composition was used 11 . 

In this article, we have calibrated the Mj-£Ca diagram using the observed values of Mj and SCa for 
stars in globular clusters, which have very old ages that are within a few Gyrs of 12.6. We have also used 
observations of the open cluster M67, which has an age of 4.0 Gyrs and a solar composition (Tautvaisienc 
ct al. 2000, and references therein) and the LMC, which has a wide range of ages. According to RGB's 
plotted in Fig. 16 for ages of 12.6 and 4.0 Gyrs, very little error in metallicity is introduced by ignoring the 
age difference between globular clusters and younger systems of the same metallicity. 

From previous investigations of the CMD of Fornax, we know that it contains stars that span the age 
range from the ages of globular clusters down to roughly 0.5 Gyrs (see Introduction). Figure 16 illustrates 
that over this age range, the Mj-SCa diagram is mostly sensitive to metallicity and only mildly to age. If, 
for example, one uses the curves in Figure 16 for 12.6 Gyrs as fiducial curves and then derives values of 
[A/H] for the stars lying along the curves for 0.25 Gyrs, without taking into account their younger ages, one 
will underestimate the true abundances of the stars by about 0.2 dex. We have tested this prediction using 
observations of the open star cluster Mil. Note that only a small fraction of the Fornax stars are likely to 
be this young and that this systematic error decreases substantially with increasing age. 

The curves in Figure 16 show that the V-I vs. SCa diagram is sensitive to age differences, particularly 
at the highest abundances where SCa is more sensitive to surface gravity. While this diagram is less useful 
for estimating metallicity than the Mj-£Ca diagram, it does provide important information when used in 
conjunction with this other diagram. Note that on the basis of the M/-SCa diagram alone, one could interpret 
stars lying on the 0.25 Gyr and [A/H]=0.0 curve at, for example M/=— 3, as old stars (age > 4 Gyrs) with 
[A/H] « —0.2. Such stars are expected to lie below the curve for [A/H]=0.0 and 4.0 Gyrs in the V-I vs 
SCa diagram. If instead they lie above this curve, the first interpretation is incorrect, and the stars are 
probably younger than 4.0 Gyrs with abundances closer to solar. 



In order to estimate ECa for the most luminous stars, we extrapolated the formulae given by Jorgensen ct al. (1992) to 3500 K, 
which is substantially beyond the coolest point on their grid, 4000 K. This may also extend the calculations beyond the point 
where the contamination of the Ca II lines by TiO becomes important. For this reason and because the observed values of 
SCa are pseudo-equivalent widths that are measured relative to nearby spectral regions that include weak absorption lines and 
therefore not the same as the true equivalent widths plotted in Figure 16, these figures cannot be used for analyzing directly the 
observed values of SCa. The curves in these figures can be used in a differential way to show the sensitivities of these diagrams 
to age and composition. 
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Fig. 1. — Position of the 117 targets on the 
plane of the sky. Smaller symbols denote tar- 
gets fainter than 1=17.5. The different symbols 
indicate: hexagons: metal-poor stars ([Fe/H]< 
— 1.3); squares: intermediate metallicities; trian- 
gles: metal-rich stars ([Fe/H]> —0.7). The star 
symbol indicates the position of the globular clus- 
ter located near the center of Fornax. 



Fig. 2. — Upper CMD of Fornax. Large dots 
indicate the measured targets and open symbols 
denote likely non-members and objects removed 
from the final analysis by the S/N criteria. The 
sequences of the globular clusters at [Fc/H] = — 
1.1 and [Fe/H]=— 0.7 (from Saviane & Rosenberg 
1999) are also indicated. 
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Fig. 5. — The positions of the red giants in Fornax (open circles) and in Mil (crosses), compared with the 
sequences of some clusters, in the (V-I)-Mj color-magnitude diagram (top left), in the SCa vs. Mj plane 
(top right), in the SCa vs. My plane (bottom left) and in the ECa vs. (V — I)o plane (bottom right). 
The lengths of the lines indicate the range covered by the red giants in the clusters. These clusters were 
used to calibrate our measurements of the Ca II triplet as a function of metallicity (see Fig. 7). 
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Fig. 3. — Sky subtraction procedure for object 
103. Top: raw spectrum. The sky emission lines 
are prominent. The three Ca II lines are hardly 
noticeable between them. Middle: the sky spec- 
trum, taken on the side of the object spectrum. 
Bottom: the final object spectrum, reduced with 
the offset method. Note the efficiency of the sky 
subtraction with this method. 
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Fig. 6. — Metallicity as a function of I magnitude 
for our sample. There is no significant increase of 
the computed metallicity with brightness, imply- 
ing that the curvature of the calibration towards 
the brightest magnitudes is limited. 
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Fig. 4. — Three representative spectra. Top: ob- 
ject with average metallicity ([Fe/H]~ —1). Mid- 
dle: object with low metallicity ([Fc/H]~ —2). 
Bottom: spectrum that is too noisy to satisfy our 
selection criteria. (Note that the average quality 
of our spectra is significantly higher than in previ- 
ous studies of Fornax such as Tolstoy et al. 2001). 
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Fig. 7. — Growth of the Ca II triplet reduced 
equivalent width with metallicity. The points in- 
dicate the calibrating clusters. The solid line is the 
cluster calibration (Equ. 1). The dotted line shows 
the modified relation for young populations with 
Mi > -3 (Equation 2), based on the LMC con- 
straint represented by the star symbol. Error bars 
are from CG (vertical) and from the propagation 
of the uncertainty on the fits of SCa vs. Mj (hor- 
izontal). Note that this figure is a 2-D projection 
of an actual relation in 4-D space (colour, magni- 
tude, metallicity, Calcium equivalent width). The 
lines in the figure are projections of the same cal- 
ibration for different magnitude ranges: Mi ~ —2 
for the clusters and Mj — -4 for the LMC. 



30 - 




[Fe/H] 



Fig. 9.— Histogram: Metallicity distribution of 
the Fornax red giant sample, based on the M/ cali- 
bration. The shaded part indicates the less secure 
part of the calibration, as explained in the text. 
Dotted line: Inferred metallicity distribution for 
the underlying population, obtained by applying 
the corrections suggested by synthetic population 
modeling (Table 3). The resulting histogram was 
smoothed with a Normal kernel of <J\FeiH\ =0.1 
dex. 
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Fig. 10. — Position of the RGB at a fixed I mag- 
nitude in the age-color plane according to the 
Yonsei-Yale isochrones. The red limit of the For- 
nax giant branch is indicated as a dashed line. 
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Fig. 11. — The age-metallicity relation of our sam- 
ple. Ages derived by comparison with stellar evo- 
lution models of Girardi et al. (2000). The hor- 
izontal bars are the 95% age confidence intervals 
for each star. The most likely age is indicated by a 
pentagon, when the confidence interval is smaller 
than 5 Gyr. The dotted line illustrates the predic- 
tion of a closed-box model with constant star for- 
mation rate and a remaining gas fraction of 0.0018 
at t = 0.5 Gyr. 
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Fig. 13. — The cumulative metallicity distribution 
for Fornax - solid line: our sample, broken line: 
underlying population - compared to a closed-box 
model and to the data of J0rgensen (2000) for the 
Galactic disc (shifted by 0.6 dex). 
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Fig. 14. — Position of the red giant branch at 
t=12.6 Gyr for different metallicities in the log Te 
vs. log g plane, with some lines of constant color 
and constant magnitude indicated. The long- 
dashed line delimits the zone where TiO bands 
stars affecting the Ca II triplet region of the spec- 
trum. 
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Fig. 15. — Red giant branch and asymptotic giant 
branch models in the log Te vs. log g for z=0.019. 
The whole red giant branch locus is indicated for 
t=12.6 Gyr, and the position of the point at Mj — 
—3 mag for different ages. 
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Fig. 12.— Left: Observed CMD for the central 6.7'x 6.7'of Fornax, obtained using FORS at the VLT, 
and published in Gallart et al. (2002). A set of lines bracketing the upper observed main sequence, and 
the Armandroff and Da Costa (1991) sequences for the globular clusters M15, M2, NGC1851 and 47 Tuc, 
have been drawn as a reference. Middle: synthetic CMD obtained assuming a constant star formation rate 
between 12 Gyr and 0.5 Gyr ago, the Kroupa et al. (1993) IMF, 25% binary stars, and the AMR shown in 
the inset, which is a fit to the AMR displayed in Figure 11. Right: same as middle panel, but modifying the 
assumed AMR to the one shown in the corresponding inset, which is approximated from the AMR derived 
by Saviane et al. (2000) from the position of the stars in the RGB (in both synthetic CMDs, a gaussian 
metallicity dispersion from the mean AMR has been introduced, with cr(z)=0.0005 at all ages). The position 
of the stars on the RGB does not match exactly in the observed and synthetic samples, and this is a known 
major drawback of current stellar evolution models when used to fit CMDs. The position of main-sequence 
stars predicted by the models, however, is much more reliable, and a good agreement is found between the 
positions of the main sequence in the observed CMD and the synthetic CMD in the middle panel, computed 
using the AMR obtained in the current spectroscopic analysis. Both main sequences are nicely centered in 
the brackets. The main sequence in the right panel, however, is displaced to the blue with respect to the 
observed main sequence, due to the lower metallicities of the stars younger that 2 Gyr. 



30 




0-2-4 -6 -2 -4 1 1.5 2 2.5 



M, M, V-I 



Fig. 16. — Ca II line strength models for red giants in the M/ vs. EW (left), My vs EW (center) and 
V-I vs EW (right) planes. For the 0.25 Gyr models, both the AGB (dashed curve) and the RGB (solid 
curve) are plotted. For other ages, only the RGB is plotted. (Note that the relations plotted cannot be used 
directly to analyse observational data - see footnote 11) 
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Table 1: Relevant data for the objects measured in the calibration globular clusters. The star designations 
follow the scheme that Rutledge et al. (1997b) used in their catalog of Ca II measurements, and the 
photometry is from the same source. 



Star ID 


V 


(B-V) 


Vobs 


EW 


a 


EW 


a 


EW 


a 


E(Ca) 


a 


E(Ca) 


a ID EC 








kms -1 


^8498 




^8592 




-^8662 




this work 




Ru97 




NGC 2298 


A14 


14.35 


1.22 


147.20 


0.79 


0.09 


2.20 


0.11 


1.77 


0.11 


3.66 


0.13 


3.46 


0.13 


A17 


14.87 


1.09 


114.60 


0.76 


0.12 


1.94 


0.13 


1.55 


0.13 


3.26 


0.17 


3.20 


0.12 


A18 


14.96 


1.00 


135.05 


0.72 


0.14 


1.86 


0.15 


1.54 


0.16 


3.14 


0.20 


3.04 


0.13 


A22 


15.30 


1.02 


127.75 


0.64 


0.15 


1.74 


0.17 


1.41 


0.17 


2.91 


0.21 






A15 


14.86 


1.12 


122.30 


0.80 


0.11 


1.89 


0.12 


1.50 


0.12 


3.19 


0.15 


3.08 


0.11 


NGC 3201 


L2207 


14.30 


1.07 


523.38 


0.87 


0.13 


2.30 


0.14 


1.88 


0.15 


3.87 


0.18 


4.00 


0.13 


L2214 


13.84 


1.03 


516.21 


0.88 


0.10 


2.17 


0.11 


1.73 


0.12 


3.65 


0.14 


3.67 


0.11 


L3204 


12.29 


1.44 


479.98 


1.26 


0.05 


3.01 


0.05 


2.35 


0.05 


5.05 


0.06 


5.08 


0.03 


L3101 


13.66 


1.16 


479.01 


1.08 


0.10 


2.67 


0.11 


2.06 


0.10 


4.44 


0.13 


4.29 


0.07 


L3102 


14.31 


1.00 


471.67 


0.99 


0.13 


2.25 


0.14 


1.78 


0.15 


3.81 


0.18 


3.80 


0.11 


L3107 


13.38 


1.14 


474.90 


0.96 


0.09 


2.44 


0.09 


1.88 


0.09 


4.05 


0.12 


4.11 


0.07 


L4215 


13.56 


1.26 


514.17 


0.95 


0.09 


2.51 


0.10 


1.99 


0.10 


4.17 


0.12 


4.32 


0.12 


L4307 


14.33 


1.11 


502.04 


0.93 


0.12 


2.33 


0.13 


1.87 


0.14 


3.92 


0.17 


3.82 


0.11 


L4318 


12.54 


1.45 


493.45 


1.17 


0.05 


2.97 


0.05 


2.35 


0.05 


4.97 


0.06 


5.03 


0.03 


L1309 


13.26 


1.03 


487.97 


0.94 


0.08 


2.48 


0.08 


1.91 


0.09 


4.09 


0.11 


4.16 


0.08 



NGC 7099 



DP24 


13.06 


1.09 


-271.84 


0.72 


0.17 


1.95 


0.19 


1.53 


0.19 


3.23 


0.24 


3.20 


0.10 


DP23 


12.20 


1.40 


-170.28 


0.86 


0.10 


2.29 


0.11 


1.84 


0.11 


3.82 


0.13 


3.64 


0.10 


DP91 


13.00 


1.06 


-122.12 


0.79 


0.18 


1.81 


0.18 


1.45 


0.18 


3.07 


0.23 


3.04 


0.10 


D31 


14.22 


0.81 


-118.27 


0.54 


0.20 


1.47 


0.22 


1.12 


0.22 


2.41 


0.27 






D57 


14.93 


0.73 


-140.11 


0.63 


0.31 


1.45 


0.32 


1.15 


0.34 


2.46 


0.41 


2.21 


0.19 


DP19 


13.14 


1.04 


-128.02 


0.79 


0.11 


1.77 


0.11 


1.41 


0.12 


3.01 


0.14 


2.91 


0.04 
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Table 2: Fornax star data. Photometry from Gallart et al. (in prep.). [Fe/H]cG is the metallicities on the 
Carretta & Gratton (1997) scale. Key to the notes: see bottom of the Table. 



ID 






alpha 






delta 




v 




crV 




R 




crR 


I 




al 




V-I 


EC a 






[Fe/H] CG 


101 


2 


40 


08 


42 


-34 


25 


33 


89 


18 


607 





005 


17 


755 





004 


17.037 





010 


1 


570 


7 


29 





24 


-0 


3: 


102 


2 


40 


06 


56 


-34 


25 


22 


57 


19 


032 





007 


18 


306 





005 


17.700 





015 


1 


332 


6 


22 





33 


-0 


6: 


103 


2 


40 


05 


07 


-34 


24 


38 


39 


18 


646 





006 


17 


917 





004 


17.299 





016 


1 


347 


4 


85 





22 


-1 


38 


104 


2 


40 


02 


76 


-34 


23 


32 


92 


18 


527 





006 


17 


645 





004 


16.868 





020 


1 


659 


6 


28 





18 


-0 


8: 


105 


2 


40 


01 


13 


-34 


25 


09 


72 


18 


208 





005 


17 


434 





003 


16.752 





020 


1 


456 


4 


78 





16 


-1 


54 


106 


2 


39 


57 


97 


-34 


23 


20 


69 


18 


425 





005 


17 


642 





003 


16.948 





027 


1 


477 


5 


32 





29 


-1 


22 


107 


2 


39 


57 


06 


-34 


23 


32 


81 




- 




- 




- 




- 


(18.287) 




- 




- 


6 


47 





42 




_ 


108 


2 


39 


54 


32 


-34 


23 


31 


10 


18 


371 





004 


17 


547 





003 


16.900 





020 


1 


471 


5 


00 





26 


-1 


40 


109 


2 


39 


53 


17 


-34 


22 


58 


75 




- 




- 




- 




- 


(17.261) 




- 




- 


3 


86 





25 


-1 


82 


110 


2 


39 


51 


33 


-34 


24 


53 


41 


18 


683 





005 


17 


876 





005 


17.235 





023 


1 


448 


5 


73 





29 


-0 


97 


111 


2 


39 


49 


78 


-34 


23 


52 


76 


18 


559 





008 


17 


846 





006 


17.256 





021 


1 


303 


4 


76 





21 


-1 


43 


112 


2 


39 


48 


58 


-34 


24 


48 


48 


19 


033 





009 


18 


263 





006 


17.883 





009 


1 


150 


5 


81 





26 


-0 


8: 


113 


2 


39 


46 


05 


-34 


25 


04 


52 


19 


040 





009 


18 


286 





004 


17.663 





017 


1 


377 


7 


14 





26 


-0 


3: 


114 


2 


39 


43 


92 


-34 


23 


59 


85 


18 


595 





005 


17 


895 
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006 
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29 


66 


18 
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006 
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798 
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279 
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18 
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41 
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785 
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17.320 
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98 
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551 
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702 
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577 
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88 
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55 
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688 
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026 
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17.417 
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1 


271 


2 


97 





23 


-2 


07 


311 


2 


39 


33 


56 


-34 


25 


49 


50 


18 


569 





006 


17 


706 





003 


16.965 





004 


1 


604 


6 


70 





21 


-0 


6: 


312 


2 


39 


35 


40 


-34 


26 


05 


90 


19 


118 





006 


18 


417 





005 


17.774 





008 


1 


344 


5 


00 





30 


-1 


17 


313 


2 


39 


32 


40 
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26 


35 


40 


19 


051 
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461 





006 


17.919 





010 


1 
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5 


87 





33 


-0 


7: 
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2 
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33 
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26 


50 


40 


18 
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18 
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1 


439 


5 


70 





23 


-0 


9: 
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2 


39 


26 


05 
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27 


06 


80 


19 
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008 


18 


434 





005 


17.805 
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1 


352 


5 


62 





30 


-0 


92 


316 


2 


39 


27 


40 


-34 


27 


27 


70 


19 


213 
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18 


460 
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17.840 





013 


1 


373 


6 


65 





32 


-0 


4: 


317 


2 


39 


21 


75 
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27 


46 


70 


18 


837 
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18 


056 
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17.387 
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1 


450 


5 


76 





24 


-0 


9: 



a[l 



Table 3: Expected metallicity bias in our sample. 



Metallicity Ratio 

range selected vs. total 



-2.0 


- -1.8 


0.80 


-1.8 


-1.6 


0.70 


-1.6 


- -1.4 


0.73 


-1.4 


- -1.2 


0.70 


-1.2 


-1.0 


0.91 


-1.0 


- -0.8 


1.00 


-0.8 


- -0.6 


1.48 


> 


-0.6 


1.15 
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Table 4: Expected age bias in our sample. 



Age 
range 


Ratio 
selected vs. total 


Proportion 
of AGB 


1-3 


1.35 


43% 


3-5 


1.35 


16% 


5-7 


0.96 


17% 


7-9 


0.93 


21% 


9-11 


0.78 


25% 


11-13 


0.73 


15% 


13-15 


0.90 


24% 



Tabic 5: Inferred age distribution and star forma- 
tion rate. 



Age 


Stars 


Corrected 


Corresponding 


interval 


in sample 


to total 


SFR 






population 


(mcan=l) 


0-2 Gyr 


35 


26 


2.4 


2-4 Gyr 


31 


23 


2.1 


> 4 Gyr 


27 


34 


0.6 
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